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Abstract 

We study the linear differential equation x = Lx in 1:1 resonance. That is, x £ R 4 and L is 4 x 4 
matrix with a semi-simple double pair of imaginary eigenvalues (i(3, —i(3, i/3, —i(3). We wish to find all 
perturbations of this linear system such that the perturbed system is stable. Since linear differential 
equations are in one to one correspondence with linear maps we translate this problem to gl(4, R). In 
this setting our aim is to determine the stability domain and the singularities of its boundary. The 
dimension of gl(4, R) is 16, therefore we first reduce the dimension as far as possible. Here we use a 
versal unfolding of L ie a transverse section of the orbit of L under the adjoint action of Gl(4, R). 
Repeating a similar procedure in the versal unfolding we are able to reduce the dimension to 4. A 
3-sphere in this 4-dimensional space contains all information about the neighborhood of L in gl(4, R). 
Considering the 3-sphere as two 3-discs glued smoothly along their common boundary we find that 
the boundary of the stability domain is contained in two right conoids, one in each 3-disc. The 
singularities of this surface are transverse self-intersections, Whitney umbrellas and an intersection 
of self-intersections where the surface has a self-tangency. A Whitney stratification of the 3-sphere 
such that the eigenvalue configurations of corresponding matrices are constant on strata allows us to 
describe the neighborhood of L and in particular identify the stability domain. 



Introduction 



1.1 Setting 

Suppose that the ordinary differential equation 

x = F M (x) (1) 

has a stationary point at x = for all /i. In this equation x € R n , /i€tt p and F^{x)-^ is a vector field on 
E™ smoothly depending on x and /i. Furthermore suppose that A : M. p — > gl(n,R) : /i i— > A(n) — DF^O) 
is the linear part of the vector field at x = 0. We are interested in the case that the phase space is 
4-dimensional and the linear part has a double semi-simple pair of complex conjugate imaginary 

eigenvalues (i/3,—i/3,i(3,—i(3), (3 ^ 0, at \i = 0. A system with such a linear part is also said to be in 
1:1- resonance. With / denoting a two by two identity matrix and (3 = 1 we have 

A(0) = L=(° I £). (2) 

Our aim is to give a description of a small neighborhood of L in gl(4, M.) and in particular find the stability 
domain and its boundary. On the latter we expect bifurcations of the stationary point of the original 
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differential equation (JTJ) . Thus this study fits in a much larger study of local bifurcations of systems 
having a stationary point with zero or imaginary eigenvalues in the linear part. 

We give a short list of generic bifurcations with increasing number of zero or imaginary eigenvalues in the 
linear part. We do not specify non-degeneracy conditions, but instead refer to the literature. In dimension 
one, only one eigenvalue can be zero and we have transcritical (TC) and saddle-node (SN) bifurcations. 
However in the presence of symmetry the pitchfork (PF) bifurcation occurs. In dimension two there can 
be two non semi-simple zero eigenvalues denoted by 2 , then the system has a Bogdanov-Takens (BT) 
bifurcation. But there can also be a pair of imaginary eigenvalues, then the system undergoes a Hopf 
bifurcation. If the two dimensional system is Hamiltonian and has two zero eigenvalues 2 , then there is 
in general a Hamiltonian saddle-node (SSN) bifurcation. 

When the dimension gets larger, results get sparser because the codimensions of the bifurcations readily 
increase. In dimension three, the simplest case is the Hopf- saddle-node (HSN) bifurcation when there 
is one zero eigenvalue and an imaginary pair, denoted 00. Other cases are systems with linearizations 
having eigenvalues 3 or 00 2 . The bifurcations for these cases are not yet well studied. The simplest 
bifurcation for general systems in dimension 4 is the Hopf-Hopf (HH) bifurcation where the linearization 
at the stationary point has two pairs of non-resonant imaginary eigenvalues 0i0i- Other results for 
4-dimensional systems are known for special cases only. For Hamiltonian systems there is a bifurcation 
at fc:/-resonance of imaginary eigenvalues 0\0%. A special place is taken by the l:±l-resonances, where 
the sign in 1:±1 designates symplectic signature of eigenvalues. At 1:— 1-resonance there is a so called 
Hamiltonian- Hopf (SH) bifurcation. But at 1: 1-resonance there is another bifurcation that has not been 
analyzed completely. The 1: 1-resonance in reversible systems is very similar to the 1:— 1-resonance in 
Hamiltonian systems, note that there is no signature for imaginary eigenvalues in reversible systems. 
However, there is in reversible systems a reversible sign for real and zero eigenvalues, see |22j . The case 
Oi has been studied, but the case 0l has not, as far as we know. We summarize the local bifurcations 
(with lowest codimension) occurring in systems up to dimension 4 in table [TJ 



dim 


eve 


bif 


codim 


comments 


references 


1 





SN 


1 










TC 


1 


u x — stationary for all //" 


[MISS] 






PF 


2 










PF 


1 


Z2-symmetric 


ma Eg 


2 


2 


BT 


2 




mi 







H 


1 




[HISS] 




2 


SSN 


1 


symplectic 


m 


3 


00 


HSN 


2 




8, 19, 17J 


4 


0102 


HH 


2 


0i : 02 irrational 


PEGS] 




2 


SH 


1 


symplectic 1:— 1-resonance 


[38] 




0102 


nn 


1 


symplectic l:2-resonance 


m 




0102 


nn 


3 


symplectic l:3-resonance 


m 




0102 


nn 


2 


symplectic k: /-resonance, kd ^ 1:2, k:l ^ 1:3 


m 




2 


nn 


1 


reversible 1 : 1-resonance 


39 




0+0+ 


nn 


nk 


symplectic 1: 1-resonance 


m 




00 


nn 


nk 


1: 1-resonance 






ol 


nn 


2 


reversible, reversible sign +1 


m 



Table 1: Low codimension bifurcations of stationary points in systems up to dimension 4- We do not list 
the non-degeneracy conditions. When such conditions are violated but higher order conditions are met, we 
generally get a higher codimension bifurcation for the same eigenvalue configuration. The abbreviations 
have the following meaning dim: dimension of phase space; eve: eigenvalue configuration ( see section 
\5.1\) : bif: name of bifurcation; codim: codimension of bifurcation; SN: saddle-node; TC: transcritical; 
PF: pitchfork; BT: Bogdanov-Takens; SSN: Hamiltonian saddle-node; HSN: Hopf- saddle-node; HH: Hopf- 
Hopf; SH: Hamiltonian Hopf; nn: no name; nk: not known. 

Here we do not even attempt to describe the bifurcation of system ([1]), whose linear part at x — and 
fi = is in 1: 1-resonance. Instead we focus on the linearized system x = A(fi)x near fi = 0. As mentioned 
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before we will describe the neighborhood of A(0) = L in gl(4,R), in particular the stability domain and 
the singularities on its boundary. Moreover it is of both theoretical and practical interest to know how 
linear Hamiltonian, reversible and equivariant subsystems appear as subspaces in gl(4, R) and how they 
intersect the stability domain and its boundary near L. This last issue will be treated in more detail in 

EH. 



The stability domain in gl(4, R) is defined as follows. 

S = {A e gl(4, R) | if A is an eigenvalue of A then Re(A) < 0} 



(3) 



Then the boundary of the stability domain <9£, is characterized by vanishing real parts of one or more 
eigenvalues of A £ gl(4, R). By the implicit function theorem, <9E as a hyper surface in gl(4, R), is smooth 
in points where the corresponding matrix has a simple zero eigenvalue or a simple complex conjugate 
pair with vanishing real part. At points where multiple eigenvalues have vanishing real parts we may 
expect singularities. A simple example being two pairs of purely imaginary eigenvalues ±i{3i and ±i/?2 
with /3i (3-2 . Then <9£ has generically a transverse self- intersection. We will only consider the stability 
domain and its boundary in a small neighborhood of L and in figure [T] we list the possible eigenvalue 
configurations. For an informal definition of eigenvalue configuration see section [2] a precise definition 
will be given in section O We use the following coding: 7 for a pair of complex conjugate eigenvalues; 
(3 for a pair of complex conjugate imaginary eigenvalues; a for a real eigenvalue; for a zero eigenvalue; 
77 for semi-simple double eigenvalues and 7 2 for double eigenvalues with nilpotent part of height 2. A 
subscript ± denotes the sign of the real part and an optional index is used to denote different eigenvalues. 
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Figure 1: Eigenvalue configurations near L. (3(3 is the eigenvalue configuration of L. On we have f3 2 , 
/?i/?2 or (3j-, on £ we only have 7-17^2 and elsewhere we have 7-7+, /?7+ or 7+17+2- 



Remark 1.1. There are many more eigenvalue configurations on dT,, including zero or real eigenvalues, 
for example ct-ia-2(3, 0a_/3 and 4 . In points of <9E where the corresponding matrix has eigenvalue 
configuration a^ia^P, dT, is smooth, for 0a_/3 it has a self-intersection and for 4 , <9£ will be more 
singular. In a study of 4 the present study of L will appear as a sub case. However, these eigenvalue 
configurations do not occur on an arbitrary small neighborhood of L. ► 

The question of singularities on the boundary of the stability domain has been taken up earlier, see for 
example a discussion of the decrement diagram in pQ. There a classification is guided by codimension 
that is by the number of parameters in a family of matrices A(/f) where ^4(0) is the central singularity. 
Also see [MJ [37] for an elaboration on this idea with examples. Here we wish to view a study of L in a 
classification guided by dimension of the phase space. This will lead to high codimensions. Indeed, in 
studying L we have to consider an eight parameter family. See section f3. 21 how we reduce such families. 
In the end it turns out that we have to study a three parameter family. In this family we do find most 
of the singularities listed in [T] for generic three parameter families. 

1.2 Motivation and main questions 

The main motivation for this study comes from a wide variety of applications where the question of 
stability of a system near l:l-resonance turns up in various forms. For example, double semi-simple 
imaginary eigenvalues are natural in the spectra of rotationally and spherically symmetrical models of 
solids and fluids [3T], ranging from car brakes [SOj, rotating shafts [3D], and computer hard discs [IT] 
to rotating elastic Earth [42], and from vortex tubes [21] to magneto-hydrodynamics 28J. Double semi- 
simple eigenvalues are also characteristic of optimal structures and are responsible for high sensitivity of 
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the latter to small imperfections [43] . Many dissipation-induced instabilities in water wave models can 
be traced back to the occurrence of double imaginary eigenvalues, see [5] and references therein, but also 
[551 [H (3j>l 1321 HS1 OH] • An early observation of friction induced instability can be found in [45 , which has 
been related to a singularity on the boundary of the stability domain by [9] . For more applications and 
references see [23] . 

The questions in the applications mentioned above in many instances boil down to questions about 
the stability domain and its boundary in gl(4, R) near a matrix with double semi-simple imaginary 
eigenvalues. The main questions we wish to address here are the following. 

1. What are the open domains in parameter space with constant eigenvalue configuration? 

2. What are the singularities on the boundaries of these domains? 

3. In particular, we address the above two questions for the stability domain. 

We will consider these questions for a general system. But in many applications such a system can be 
considered as small, dissipative perturbations of a Hamiltonian system in 1 : ±l-resonance. However, also 
reversibility and equivariance with respect to a circle group play a prominent role. In all our constructions 
to study the neighborhood of a matrix with double semi-simple imaginary eigenvalues, we take care to 
carry them out in such a way that Hamiltonian, reversible and equivariant subsystems can be recognized 
easily. 

1.3 Organization 

In section[5]we give an overview of the results, that is a description of the domains with constant eigenvalue 
configuration near L in gl(4,R). In particular we present the stability domain and its singularities. The 
methods we use will be outlined in section [3J We give a short overview of the centralizer unfolding of a 
matrix A in gl(n,R), a family describing a neighborhood of A. We also present a method to reduce the 
number of parameters in this family leading to a reduced centralizer unfolding with the same properties 
but easier to analyze. This will be applied in section Q] to the matrix L introduced in section ITTTI The 
resulting reduced centralizer unfolding of L is further analyzed in section [5] where we characterize the 
stability domain and its singularities. Finally in section [5] we give a description of a small neighborhood 
of L in gl(4,R). 

2 Statement of results 

Our aim is to describe a neighborhood of a matrix with a semi-simple double pair of imaginary eigenvalues 
in gl(4, R). The dimension of this space is rather large, namely 16, therefore we first apply several reduc- 
tions to simplify the analysis. Two main ingredients of this reduction are transversality and equivalence 
classes. In the following we use some results for smooth group actions. As a general reference see [5] and 
[13] . for application to dynamical systems see [1]. 

We use the fact that a linear differential equation transforms like a linear map under a change of co- 
ordinates. That is the map A E gl(4,R) in x = Ax transforms like A i— ► g~ x Ag under the coordinate 
transformation g 6 G1(4,R). Thus, in this respect G1(4,R) acts on gl(4,R) by similarity transforma- 
tions, the action is called the adjoint action of G1(4,R) on gl(4,R). The orbits of the adjoint action are 
smooth manifolds, but moreover they form equivalence classes. 

To reduce a neighborhood of L we consider a complement of the equivalence class of L. This runs 
as follows. In view of the previous, the equivalence class of L has a tangent space. We now take a 
complement of the tangent space in L, ie locally a transverse section of the Gl(4, R)-orbit of L, see figure 
[2j In the present context this is also called a versal unfolding of L. Then for each A e gl(4, R) in a small 
neighborhood of L, the Gl(4, R)-orbit of A transversally intersects this complement. Thus each element 
A in a neighborhood of L is equivalent to an element in the transverse section of the Gl(4, R)-orbit of L. 

The results below are formulated for a specific versal unfolding of the matrix L as defined in equation ^ . 
However we only list properties of L and a neighborhood of L in gl(4, R) that are invariant under smooth 
changes of coordinates in both the phase space R 4 and gl(4, R). This means that in the neighborhood of 



Stabdom A..., Hvn, Krv..., 2009/nov/6 



5 



any other matrix in gl(4,R) with eigenvalues (if3,—if3,i{3,—i/3), (3 ^ 0, there is a stability domain with 
singularities on its boundary diffeomorphic to the one described below. 




Figure 2: Schematic picture o/gl(4, R) in a neighborhood of L. The centralizer unfolding £ is an orthog- 
onal complement (a particular transverse section) of the tangent space of Orb G \(4 y R) (L) at L. For A e £, 
Orbc\ L (4n) (A) is the G1l(4, M.)-orbit of A in £. The reduced centralizer unfolding C r is a transverse 
section of the Gli(4, R) -orbits in £. 

There are many choices for a versal unfolding. One with particular nice properties is the centralizer 
unfolding, see section I37T1 In the following result we use a basis in gl(4, R), for details we refer to section 

ED 

Lemma 2.1. Let (Mi, . . . , Mg, Pi, . . . , Pg) be a basis for gl(4, R), such that (Pi, . . . , Pg) spans the tangent 
space of the Gl(4, WL)- orbit of L and (Mi, . . . , Mg) spans its orthogonal complement. Then the centralizer 
unfolding of L is given by £ : R 8 — > gl(4, R) with 

8 

C((j,) = L + Y,ViMi. 

i=l 

Remark 2.1. In choosing a basis for gl(4,R) we have taken into account that we wish to recognize 
Hamiltonian, reversible and equivariant subsystems. ► 

By taking a versal unfolding of L we have reduced the dimension of the original problem from sixteen 
to eight. The latter is the codimension of the tangent space of the Gl(4, R)-orbit of L. One of the 
nice properties of the centralizer unfolding is that there is an easy characterization of the subgroup of 
G1(4,R) preserving £, namely G1l(4,R) the group of all transformations g commuting with L. The 
Gli(4, R)-orbits in £ form equivalence classes in £ which allows us to further reduce the dimension of 
the problem by again taking a transverse section to these orbits. This results in a reduced centralizer 
unfolding, reducing the dimension of problem by three. For details we refer to sections 13.21 and 14.21 

Lemma 2.2. A reduced centralizer unfolding of L is given by £ r : R 5 — > gl(4,R) with 
C r {v) — L + v x Mi + v 2 M A + v 3 M 6 + z/ 4 M 8 + u 5 M 5 . 

See figure [5] for a schematic picture of the relation between £ and £ r in gl(4, R). 

Remark 2.2. There are many possible choices for a reduced centralizer unfolding. This particular choice 
enables us to easily recognize Hamiltonian, reversible and equivariant subsystems later on. ► 

Instead of describing the stability domain in the full 16-dimensional neighborhood of L we do this for 
the reduced unfolding C r of L, which is still 5-dimensional. The boundary of the stability domain is 
determined by the property that L r (v) has at least one pair of purely imaginary eigenvalues. Note that 
there are no zero eigenvalues in an arbitrary small neighborhood of L. This leads to a condition for v 
and we find the following, see section l5Tl 
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Lemma 2.3. The boundary of the stability domain of C r is contained in the critical set C = \v G 
M 5 | F{v) = 0} where F{u) = {v\ - v$){v\ + v\) + v\v\. 

Note that F does not depend on v§. This is not surprising when we mention that M5 = L. Thus the 
problem of describing the stability domain has been reduced to 4 dimensions. One further reduction is 
possible when we use the fact that F is a homogeneous polynomial. This implies that the boundary of 
the stability domain transversally intersects the 3-sphere \ \v\ \ = r. Therefore we restrict the description 
of the stability domain to this 3-sphere. This only holds for small r, but after an appropriate scaling we 
may set r = 1. Then we have the following results, see section fSTTT 

Lemma 2.4. The boundary of the stability domain of the reduced centralizer unfolding C r transversally 
intersects the 3-sphere in the critical surface S = {v S M 4 | F(y) — and \\v\ \ = 1}. The critical set C 
is a (straight) cone over the critical surface S. 

With help of a parameterization of the critical surface S we obtain a Whitney stratification of S and thus 
of the 3-sphere = 1. On the strata the eigenvalue configuration is constant, so that we are able to 
identify the stability domain. For a precise definition of eigenvalue configuration see section 15.11 here we 
give an informal definition. The relevant eigenvalue configurations are shown in figure [TJ 

Definition 2.1. An eigenvalue configuration is informally defined as an equivalence class of sets of eigen- 
values such that in any two sets from the same class there are corresponding eigenvalues with negative, 
zero or positive real parts. 

We consider the 3-sphere as two 3-discs, labelled D + and D~ , glued smoothly along their boundaries, 
which are 2-spheres. Then the critical surface also has two parts S + and S~ , one in each 3-disc. We 
do this in such a way that the critical surface transversely intersects the boundary of the 3-disc. The 



parameterization of S ,± 
section 1 5.1 




[0,2tt] -» S ± : (s,t) 



: (s, i). For details we refer to 




t 



Figure 3: Domains of the map ± . On the left the domain for S + , on the right the domain for S~ . 
The images of the domains are glued smoothly along the dashed lines according to the arrows. Lines 
with three and four arrows correspond to the great circles drawn in figure [51 The 2- dimensional open 
parts are denoted Si, ... , Si, the 1-dimensional open parts L\, . . . , Lq. There are self-intersections along 
L±, . . -,Lq. Singular points are P$ and Pq where lines of self-intersection intersect. The points P\, . . . , P4 
correspond to Whitney umbrellas. On the thick lines (f> is two-to-one, therefore lines and points with the 
same name are identified in the image. 



The domains of the maps (f> are shown in figure [3] In the domains we have points Pi , . . . , Pq on lines 
Li, . . . , Lq. On these lines the maps ± are two-to-one. To avoid heavy notation we denote the images 
under ± of the various parts in the domains by the same names. The boundaries of the domains (dashed 
lines in figure [2]) are glued by the maps (jy^ , according to the arrows in figure [3] From this we see that for 
example the open 2-dimensional part S2 is bounded by P5 U L5 U P$ U Lq whose image is a (great) circle 
in the 3-sphere. Thus the image of S2 is a topological 2-disc. The same conclusion holds for Si, S3 and 
S4. Furthermore Si and S4 are glued along the circle P5 U L 5 U P& U L 6 forming a topological 2-sphcrc 
which encloses a topological 3-disc. The latter we call Vi. The 3-discs V2, V3 and V4 are constructed in 
the same way, also see figure [4] 
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This leads to a Whitney stratification of the 3-sphere. If we consider points as 0-discs, all strata are 
(topological) discs. The strata are V\, . . . , V4, Si, ■ ■ . , S4, Li, ... ,Lq, Pi, . . . , P$. In figure [5] the strata on 
D + are shown. A similar figure for D~ can be given, see figure [7] The three dimensional strata are not 
explicitly drawn in the figures. We have the following result on the stratification of the 3-sphere and the 
critical surface S. 

Theorem 2.5. The collection {Vi, . . . , Vi, Si, . . . , S4, L±, . . . , Lq, Pi, . . . , Pq} is a Whitney stratification of 
the 3-sphere. The incidence diagram is given in figure^ The sub collection {Si, . . . , S4, Li, . . . , Lq, Pi, . . . , Pq} 
is a Whitney stratification of the critical surface S . The eigenvalue configuration on the strata is given 
in tabled In particular we find that V3 is the stability domain on the 3-sphere. 




0-dim 



1-dim 



2-dim 



3-dim 



Figure 4: Incidence diagram of the stratification of the 3-sphere determined by the critical surface S. Only 
incidences of n- and n+1- dimensional strata are shown. The solid lines form the simplified incidence 
diagram used to determine the global structure of S. All strata are discs. 

Remark 2.3. The dashed lines in the incidence diagram connect points and lines to the surfaces Si 
that are not important for the global structure of the critical surface. The reason is that for example 
Pi U Li U P5 is a line segment on Si and S3 and Pi is only connected to L\. This means that shrinking 
Li, so that Pi goes to P5, does not change global connections. ► 

Remark 2.4. Since a bifurcation set is in general a semi-algebraic set it admits a Whitney stratification, 
see for example [37] ■ This is of great help in organizing the parameter space. On each stratum, systems 
are similar in some sense. To be more precise in our example on each stratum, systems have the same 
eigenvalue configuration. For an example where the organization of the parameter space is far more 
intricate and how the Whitney stratification provides structure, see [Jj. ► 
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Table 2: Eigenvalue configurations on the strata Vi, . . . , V4, Si, . . . , Si, L\, . . . , Lq, Pi, . . . ,Pg near L. 
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Figure 5: Singularities of the critical surface S. Left a Whitney umbrella with standard form x 2 = yz 2 . 
Right an intersection of self-intersections with standard form z 2 — x 2 y 2 (or equivalently z 2 = xyz). 

So far for the global structure of the critical surface. On the critical surface S we find several singular 
points and curves. In figure [5] we show the singularities at the points P\, . . . , P4 and P5 and Pq. 

Theorem 2.6. On the lines L\, . . . , Lq the critical surface S has a transverse self-intersection. In the 
points Pi, . . . ,Pq the self-intersection is non-transverse. In each of the points Px,...,Pa the critical sur- 
face S has a singularity called Whitney umbrella. But in the points P5 and P§ curves of self-intersection 
transversely intersect and the critical surface S has a self-intersection with coinciding tangent planes of 
saddle type. 

The stability domain of the reduced centralizer unfolding is V3 . Let us summarize the singularities on the 
boundary of V3 in connection with the classification in pQ. On the lines L\, . . . ,Lg we have transverse 
self-intersections with standard form xy = 0. In each of the points Pi, . . . , P4 we have a Whitney umbrella 
with standard form x 2 — yz 2 and in the points P5 and Pq we have an intersection of self-intersections 
with standard form z 2 = x 2 y 2 . 

We conclude with two pictures showing both global and local aspects of the critical surface S. Figure 
H2 shows the restriction S + of the critical surface to one of the 3-discs, namely to D + . In fact it is a 
projection onto one of the coordinate hyper planes (z^ = 0). In figure [7] we show a diagram of the global 
connection of the conoids. Since the critical surface S is on the 3-sphere we can only give a schematic 
picture. 

Remark 2.5. It is worth noting that as the image of <fi + , S + (and also S~ as the image of <j)~) is a ruled 
surface. The rules are parallel to the v\, ^4-plane and all of them pass through a curve, the ^2-axis. Such 
a surface is called a Catalan surface. In our case the curve is a straight line, then the surface is called a 
right conoid. In particular S + is known as Pliicker's conoid for n = 1, where n is the index of the family 
(s, t) 1— > (s cos t, cos nt, s sin t) of parameterizations of Plucker conoids, see [3], [T5] . ► 



3 Methods 

3.1 The centralizer unfolding 

Given a map L € gl(V') where V is a real vector space. We are interested in all maps in an open 
neighborhood of L in gl(V). Some of those will be elements of the Gl(y)-orbit 

Orb G1(v) (i) = {g- x Lg \ g e Gl(V)} 

of L and are thus equivalent to/with L. A systematic way to explore the neighborhood of L in gl(V) is 
given by Arnol'd pQ. Here we follow the description given in [53]. At several places we will tacitly use 
the fact the Gl(V)-orbit of L e gl(V) is a smooth subset of gl(V), see [51 [T3]. 

Definition 3.1. A smooth map L : W — > gl(V^) : /1 1— > £(a0 with C(0) — L is called an unfolding or a 
deformation of L. If £ is transverse to the G\{V)-orbit through L at L, then it is said to be versal. 
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Figure 6: On the left the surface S + , a right conoid, consisting of the open 2- dimensional strata Si, . . . , S4 
is shown. On the right we separately draw the lines of self-intersection L\, L2, L§ and Lq and the singular 
points Pi, P2 and Pg. At points P\ and P2 we find Whitney umbrellas and at P5 we find an intersection 
of lines of self-intersection. 




Figure 7: Diagram of the global connection of the right conoids S + and S~ . On the left the two right 
conoids are shown together with their connection along a circle of self-intersections. On the right we show 
the singular curves and points. The circle P5 U L5 U Pq U Lq is a great circle on the 3-sphere and also the 
circle containing the segments Pi U Li U P5 U L2 U P2 and P3 U L3 U Pq U L4 U P4 is a great circle on the 
3- sphere. In the points Pi, P2, P3 and P4 we have Whitney umbrellas and in the points P5 and Pq we 
have intersections of self-intersections.. 
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We are especially interested in versal unfoldings having a minimum number of parameters. The main 
source of these concepts and ideas is Arnol'd pQ, which we review in the following. 

Definition 3.2. Two unfoldings A(n) and B(p) of L are called equivalent if they are similar as families 
of linear maps. This means that there is a smooth family of transformations g(jJ,) € G1(V) such that 
g(n)A(n)g(n)~ 1 = B(p) for all /i S W. An unfolding C of L is called miniversal if (a) for every other 
unfolding A : W — > gl( V") of L there exists a smooth map x '■ ^ 9 ~~ * ^ p such that A is equivalent to Cox, 
and (b) C has the minimal number of parameters possible for unfoldings with this property. 

The number of parameters for a miniversal unfolding is equal to the codimension of the Gl(V)-orbit 
through L and so is called the codimension of L. Arnol'd [1] showed that miniversal unfoldings can 
be obtained by taking orthogonal complements to tangent spaces of Gl(V)-orbits. Such unfoldings are 
called centralizer unfoldings. To define these we first need an inner product on gl(V). The proof of the 
following lemma is rather straightforward. 

Lemma 3.1. The bilinear form (A,B) = Trace(A*B). for A, B 6 gl(V), is an inner product. 
With this inner product we have the next result. 

Proposition 3.2. The subset {L + M* | [L, M] = 0} of gl(V) is a miniversal unfolding of L. 

Proof. A short computation shows that M is an element of the orthogonal complement of the tangent 
space of Orb G i(v) (L) if [L, M*\ = 0, where [A, B] = AB - BA. U 

Let m = {M £ gl(V) | [L, M] = 0} then it is easily seen that m is a Lie-algebra. We now define the 
centralizer unfolding as follows. 

Definition 3.3. The centralizer unfolding of L is {L + M* Mem}. 

By applying the adjoint action of G1(V) on gl(V) to this unfolding we obtain an unfolding at any other 
point on the Gl(V)-orbit through L. Transversality and miniversality are preserved by this transforma- 
tion, but orthogonality will usually be lost. 

The centralizer unfolding has nice properties. Let G\l{V) = {g £ G1(V) | gL = Lg}, then the adjoint 
action of G1l(V) preserves m. Consequently there is a Lie-subgroup M. C G\l{V) such that the adjoint 
action of M preserves the centralizer unfolding. This property allows us to find equivalence classes in C 

Remark 3.1. The construction we just described can also be applied to any Lie sub-algebra g C gl(V) 
with the associated Lie subgroup G C G1(V) as transformation group. We can further generalize this to 
subsets g of gl(F) that are not necessarily Lie algebras but still preserved by the adjoint action of some 
Lie subgroup G of G1(V), see [24]. ► 

Remark 3.2. In the example we are considering in this article we have L* = —L. In a slightly more 
general situation when L* = ±L we have M. = G1l(V). ► 



3.2 The reduced centralizer unfolding 

The centralizer unfolding of L is the orthogonal complement of the tangent space of Orbci(v)(L). In 
fact, as we already noted in section I3TT1 it is a Lie sub-algebra m of gl(V). Then there is a Lie subgroup, 
which we denote by Ai, of G1(V) which acts on this complement. The Lie algebra of M. is m. This 
means that m is invariant under the Ad^-action ij!):Aixra->ni:(j,A)^ g~ 1 Ag. Thus it is possible 
to find equivalence classes in the unfolding of L. The orbit space of this group will have a lower dimension 
than the original unfolding, thus simplifying the analysis. 

However, there are two reasons not to take the quotient with the full group. The first is that it might 
not be compact so that the resulting quotient space is no longer Hausdorff. Therefore we should at least 
restrict to a compact subgroup (preferably the largest). The second reason is that even when restricting 
to a compact subgroup, the quotient space may have singularities if the action <j> is non-free. Since our 
goal is to find the singularities of the boundary of the stability domain in parameter space, we want to 
avoid introducing singularities by taking a quotient. A general theorem for compact group actions tells 
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us that the orbit space of a smooth, proper and free action (p is again a smooth manifold, see [51 [13]. In 
our case we call the quotient space m/ ~ <f> the reduced unfolding of L. The details very much depend on 
L. Therefore we refer to section [32] for the group M, its maximal compact subgroup and the quotient 
space m/ ~ <p for our particular choice of L. 

3.3 Characterization of the stability domain 

The stability domain E near L is an open set in gl(4,R), which is most easily characterized by its 
boundary <9£. However, we do not directly characterize 9£, but instead determine the set of A £ gl(4, R) 
near L where A has at least one pair of imaginary eigenvalues. We describe our method for n dimensions, 
but slightly specialized for our situation. 

First we define a map ip from sets of roots to real polynomials exactly having these roots. The map ip 
parameterizes surfaces in the space of polynomials. Then we characterize the surface where polynomials 
have at least one pair of imaginary roots as the zero set of a function /. Next we define a map <p, mapping 
elements of the unfolding of L to eigenvalue polynomials. Requiring that the eigenvalue polynomials lie 
on the surface parameterized by ip, we obtain a surface characterized as the zero set of fo p. This surface 
contains the boundary of the stability domain, see figure [5] 

ip A 
W 1 - 1 ► R n+1 ■< MP 

/ 

T 

R 

Figure 8: ip maps sets of roots with at least one imaginary pair to n-th degree polynomials, <b maps 
elements of an unfolding to n-th degree eigenvalue polynomials. The zero set of f consists of polynomials 
with at least one imaginary root pair. 

We identify the sets of roots containing at least one imaginary pair with R n_1 . The space of nth- 
degree real polynomials can be identified with R n+1 , namely a = (a n , ... ,ao) £ R" +1 is in one to one 
correspondence with p(A) = a n X n + a„_iA ,l_1 + • ■ • + aiA + oq. Furthermore we identify the p-parameter 
homogeneous (reduced centralizer) unfolding Ti{v) of L with R p , see section |4~T1 for details. 

Now we define the map ip : R™ _1 — ► R™ +1 by ib(cr) = a. This map is quasi-homogeneous because the 
coefficients a are homogeneous polynomials of the roots with different degrees. In particular if ip(a) = a, 
we have for all t £ R>o 

ip(ta) = (t n a ,t n ~ Oi, . . . ,ta„_i,a„). 

We will usually set a n = 1, then the map ip parameterizes a hyper-surface in R™. The latter represents 
polynomials with at least one imaginary pair as roots. Using the Buchberger algorithm [13] we can 
eliminate a from ip(o~) = a and obtain an implicit equation f(a) — where / : R Tl+1 — > R is a homogeneous 
polynomial. If ip{o~) is an element of the hyper-surface, then ip(to~) is also an element of the hyper-surface 
for all t £ R>o- Therefore / is also quasi- homogeneous and we have for all t £ R>o 

f(t n a , t^ai, . . . , tOn-i, On) = t k f(a), 

for some integer k. 

Let <p : MP ^ R ra+1 be the map that maps unfolding parameters v to the coefficients a £ R ,l+1 of the 
eigenvalue polynomial of Ti.(i>). Then p is quasi-homogeneous and in particular if <p{v) = a, we have for 
all t £ R >0 

<f>{tv) = {t n a , t n ~ 1 a 1 , . . . , ta n -i, a n ). 

Note that since 4>(v) are the coefficients of an eigenvalue polynomial, we always have a n = 1. 

In order that Tl(v) has at least one pair of imaginary eigenvalues, the coefficients a of its eigenvalue 
polynomial must lie on the hyper-surface defined by f(a) = 0. Thus wc find that the parameters v must 
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satisfy f{4>(v)) — 0. Therefore we now define F : M p — > R as 
F = focf>. 

The zero set of F determines a hyper-surface in the homogeneous (reduced centralizer) unfolding of L, 
which contains the stability boundary. From the homogeneity properties of ijj, f and 4> it follows that F 
is homogeneous of degree k, with the same k as above, for all t £ R>o 

F{tv) = f{cf>{tv))=t k F{v). 

Remark 3.3. On the boundary of the stability domain, a defined by <f>(v) = a must satisfy f(a) = and 
a n = 1. For our purposes it is not enough to look at the singularities of the hyper surface f(a) = because 
in the space of coefficients of eigenvalue polynomials we only have algebraic information of the eigenvalues. 
In the space of unfolding parameters we also have geometric information about eigenvalues. Thus the 
hyper surface defined by F{v) = has possibly more singularities than the hyper surface /(a) = 0. For 
example, in the latter we will not sec the difference between semi-simple and non-semi-simple eigenvalues. 
► 

Remark 3.4. In view of the previous remark it would be interesting to know the fibers of the map cj>. ► 



4 Unfolding of L 

4.1 The centralizer unfolding of L 

Here we apply the results of section l3~T1 to the map L S gl(4,R) from section [TJ which has a double pair 
of eigenvalues ±i. From the zero commutator criterion in proposition 13 . 21 the centralizer unfolding of L 
is readily seen to be 

C(n)=(° I ) + ( A ^ 
W {-I 0j + \-B(ri A(p)J- 

A(fi) and £?(/i) are any 2x2 matrices. Thus the centralizer unfolding has eight parameters and conse- 
quently every other miniversal unfolding of L must contain eight parameters. To further characterize the 
unfolding we now choose a basis for m = {M £ gl(4,R) | [L, M] = 0} and gl(4, R), regarded as a linear 
space. First we introduce 

?)■ *=g °)< t =£ s> '-(a ;)< < 4) 

as a basis for the space of 2 x 2 matrices, with help of these we define the following 16 matrices 



Mo.)< 



(5) 



i oy ' v -i? 0/' ' v-T 0/' ° vo J 

and 

Pl = (o -V) ' Pa = (o -°r) ' Ps = (o -°r) ' Pi = (j o 

Ps = (? o) ' Pe = (r o) ' Py = (t o) ' Ps = (o -J y 

Note that Mi is equal to the identity and Me, — L. The following properties are easily checked. 
Lemma 4.1. The set (M 1; . . . , M 8 , Pi, ... , P 8 ) «s a basis for gl(4,R). 



(6) 
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Lemma 4.2. The set (Mi, . . . , Mg) zs a basis for m. 

With these two lemmas we immediately get the main result of this section about the centralizer unfolding 
of L. 

Proposition 4.3. C(fi) = L + Ylj—i f^jMj, with parameters fii,...,fis € K and matrices Mi as in 
equation ([5]) 7 is a centralizer unfolding of L. The codimension of the unfolding is 8. 

The centralizer unfolding of L can be regarded as a linear variety in gl(4, K), the sum of the fixed vector 
L and the linear space m. It will turn out to be useful to have the following notion. 

Definition 4.1. We will call 7i(/i) = X)i=i fcMi, with parameters . . . , [ig € M and matrices Mi as 
in equation ([5]), the homogeneous unfolding of L. 

Remark 4.1. The matrices Mi in the proposition are not unique, since they reflect the choice of a basis in 
the vector space m. However, this particular choice will turn out to be convenient in later computations. 
► 

Remark 4.2. Note that Mb = L. When we introduce new parameters by setting /15 = vq + and 
Vi = fii for i =/= 5 we obtain C(v) = (1 + vq)T + Ti(v). Thus C is in fact the unfolding of the family of 
matrices (1 + vq)L provided that 1 + vq 7^ 0. From this we may already infer that the unfolding parameter 
/L/,5 will be relatively unimportant, also see section 15. II ► 

Remark 4.3. In connection with the previous remark, we may also set = /X5 + 1 and Vi = \Xi for 
i 7^ 5, then C(p) = Tiif). The homogeneous unfolding has the advantage that TL(tu) = tJiiv) for (el. 
Therefore, since A and tA are equivalent eigenvalue configurations for t £ K>o, and 7i(tv) have 

equal eigenvalue configurations. However TL^) is an unfolding of L only if 1/5 ^ 0. ► 

4.2 The reduced centralizer unfolding of L 

The aim of this section is to reduce the number of parameters in the centralizer unfolding of L. By this we 
mean that each element of the centralizer unfolding is equivalent to an clement of the so called reduced 
unfolding with less parameters. Here equivalence is determined by parameter dependent coordinate 
changes. As explained in section 13.21 this is equivalent to taking the quotient of m with respect to the 
adjoint action of the Lie group M. corresponding to m. In our specific example L* = — L and therefore 
M. = G1l(4,]R). Because we wish to find the singularities of the boundary of the stability domain of L 
we take care to avoid quotient singularities. Thus we will take a suitable subgroup of M. such that the 
quotient space is again a smooth manifold. 

Since we reconstruct M. from its Lie algebra m we have enough information about M. to find all elements 
with compact orbits. A table of commutators of m allows us to identify the largest appropriate subgroup 
of A4. The result is given in the next proposition. 

Proposition 4.4. A reduced centralizer unfolding is given by C r : R 5 — > gl(4,]R) and 

C r (v) = L + viMi + u 2 M 4 + iy 3 M 6 + u A M s + u 5 M 5 . (7) 

Remark 4.4. The reduced unfolding is not unique, since it represents an equivalence class. We may 
also take L + viM\ + V2M4 + i/^Mr + i^Mg + v^M^ for example. The choice made in the proposition is 
convenient when considering symplectic, reversible and equivariant subsystems, see |23j . ► 

Proof. In our setting we consider the adjoint action of M. on m 

g ■ A = g~ l Ag, 
for g € M. and A S m. 

Largest subgroup of M. with compact orbits. Every element A £ m can be written as 

8 

A = ^a i M l , 
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for certain on € R. Furthermore every element of M. (in the component continuously connected to the 
identity) can be written as exp(^* =1 ij-Mj), with t\ S R. Let us therefore compute exp(iMj) for each i, 
then we have exp(tMi) = Ci(t)id + Si(t)Mi, where 

a) ci(t) = exp(i) and s 4 (t) = 0, 

b) Ci(t) = cosh(t) and Sj(t) = sinh(i) for i e {2,3,4}, 

c) Ci(t) = cos(i) and Sj(t) = sin(i) for i S {5, 6, 7, 8}. 

Thus M5, Mg, M7 and M§ generate compact orbits. However, since M5 = L the adjoint action of 
exp(tMs) on m is trivial and the same holds for M\ = id. With help of the table of commutators, see 
table Owe identify the largest subgroup of M. with compact orbits. 





M 2 


M 3 


M 4 


M 6 


M 7 


M 8 


M 2 





2M S 


2M 7 





2M 4 


2M 3 


M 3 


-2M$ 





-2M 6 


-2Af 4 





~2M 2 


M 4 


-2M 7 


2M 6 





2M 3 


-2M 2 





M 6 





2M 4 


-2M 3 





-2M S 


2M 7 


M 7 


-2M 4 





2M 2 


2M S 





~2M 6 


M s 


-2M 3 


2M 2 





-2M 7 


2M 6 






Table 3: Commutators [A, B] in m, A vertical and B horizontal. Note that M 4 = id and A/5 = L 
commute with all elements of m. 





M 2 


M 3 


M 4 


M 6 


M 7 


M s 


AdtM 6 


M 2 


cM 3 - sM 4 


cM 4 + sM 3 


M 6 


cM 7 + sM 8 


cM & - sM 7 


Ad t M 7 


cM 2 + sM 4 


M 3 


cM 4 - sM 2 


cM e - sM s 


M 7 


cM 8 + sM 6 


AdtM 8 


cM 2 + sM 3 


cM 3 - sM 2 


M 4 


cM 6 + sM 7 


cM 7 - sM 6 


Ms 



Table 4: Adjoint action of Lie group Mers on m - Notation: c = cos(2t) and s — sin(2i). 

Note that there are several Lie sub-algebras in m, for example m678 = {M = aM^ + bMr + cMg | a,b,c € 
R} which is similar to so(3). We denote the corresponding Lie group by Mars- Clearly M^rs is the 
maximal Lie subgroup of M. with compact orbits acting non-trivially on m. 

Action of ^678 on m - Let us compute the adjoint action of this subgroup in tabled] by using 
AdtAi(A) = e~ tM Ae tM 

with the generators AdtM 6 , Ad r M 7 and AdtM 8 of Mers- Then the adjoint action of Mq7s on m is an 
SO(3)-action. Using the fact that every A £ m can uniquely be written as A — 53i=i a i^i we now 
obtain an SO(3)-action on R 8 . Again from tabled we have 

8 8 

Ad tMk (A) = e- tM *Ae tM * =Y l ^- tM "M i e tM " =^^((,0)^, 

i=l i=l 

where IZk is defined with help of the rotations Ri, R 2 and R 3 in R 3 around the x±, x 2 and x 3 axes. 
Furthermore when we write a — (ai, x, a§, y) G R x R 3 x R x R 3 and x = (a 2 , a 3 , a 4 ), y — (q.q, a?, as) 
then 

TZ 6 (t,a) = (a 1 ,R 1 (-2t)x,a 5 ,R 1 (2t) y), 
TZr(t,a) = (ai,R 2 (-2t)x,a 5 ,R 2 (2t) y), 
TZs(t,a) = {ai,R 3 ( 2t) x, a 5 , R 3 (2t) y). 
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The action on R 8 is generated by TZq, IZ7 and TZs- As an SO(s)-action it preserves ||a|| so we may 
restrict to the 7-sphere ||a|| = 1. Now we see that the action of Aie7s acts trivially in the ct\ and a.5 
directions. Thus for all fixed values of ct\ and as we have a non-trivial action on the 5-sphere defined by 
{(x,y) G M 3 x R 3 I ||a;|| 2 + \\y\\ 2 = 1}. Considered as an action on M 6 = R 3 x M 3 , K acts as the sum of 
two isomorphic standard representations of SO(3) on R 3 . 

Orbits and orbit types of the A^678 _ac tion. The orbits of the ./Vf678-action on S 5 are two or three 
dimensional. Points with three dimensional orbit have trivial isotropy group whereas points with two 
dimensional orbit have isotropy group SO(2). To see this we consider the tangent space of the A^s-orbit 
of (x, y). It is spanned by the three vectors (— L±x, L±y), (—L2X, Liy) and (L32;, L^y). Where L\, L2 and 
L3 are the standard generators of so(3). These vectors are linearly independent provided that y ^ Xx for 
A G R, then the tangent space is three dimensional. If (x,y) = (az,(3z) with a 2 + /3 2 = 1 and \ \z\ \ = 1 
then only two vectors are linearly independent, in this case the tangent space is two dimensional. 

The action of 1Z is smooth and proper but it is non-free. Therefore we can not use the general theorem 
about smooth orbit spaces referred to in section 13.21 The points on the 5-sphere with a non-trivial 
isotropy group are {(ax, fix) e I 3 x I 3 ||x|| 2 — l,a 2 + /3 2 = 1} and the isotropy group is SO(2). This 
is the only non-trivial isotropy group, so there are two orbit types SO(3) and SO(3)/SO(2). Then it is 
a result of |H] and in a somewhat broader context [35] that the orbit space m/ ~ <f> is a 2-disc. Here <f> 
is the action of M. 678 on S 5 . 

Computing the quotient. Since the action 1Z is explicitly given, we can perform explicit computations. 
For every (x, y) G S 5 , that is | \x\ \ 2 + | \y\ \ 2 — 1, we can find a rotation in SO(3) taking (x, y) into (£363, y). 
Here e% is the third standard basis vector in R 3 , £3 = ||x|| and y is the rotated vector y. If we now apply 
TZ%, leaving e% invariant, we can take (£363, y) into for example (£363, 771 e\ + 77363) with £| + nf + rj 2 = 1 
and 771 > 0. Since £3 > and 771 > this defines a 2-disc, although with a non-smooth boundary. We 
do not want a non-smooth quotient space nor do we want a boundary, therefore we take the quotient 
of S 5 with an SO(3)/(Z2 x Z2)-action so that the orbit space becomes a smooth 2-sphcre, defined by 
{(x,y) g R 3 x R 3 I £f + T)\ + v 2 3 = l,^ = 6 = % = 0}. ■ 

5 Stability domain of the reduced unfolding of L 

The main goal of this section is to find the stability domain E in parameter space of the reduced unfolding 
of L. In particular we are interested in the boundary of S. Furthermore we wish to locate several 
eigenvalue configurations which are relevant for systems with additional structure. 

5.1 Characterization of the boundary of the stability domain 

The stability domain and in particular its boundary is defined in terms of eigenvalue configurations rather 
than numerical eigenvalues. On the boundary the corresponding maps have at least one eigenvalue 
zero or a single pair of purely imaginary eigenvalues. The map L has a double pair of non-zero semi 
simple imaginary eigenvalues, therefore not all eigenvalue configurations occur on an arbitrary small 
neighborhood of L. Indeed by continuity there can be no zero eigenvalues arbitrary close to L. Thus only 
the eigenvalue configurations (3j3, (3 2 , and f3j- are relevant for the boundary, see figure [1] 

Let us now be more precise about what we mean by eigenvalue configuration. It is an equivalence class 
similar to that in |34j but defined in a different way. 

Definition 5.1. Two collections of eigenvalues A = (Ai,--- ,A„) and A = (Si,-- - ,S n ) belong to the 
same eigenvalue configuration if a permutation tt of {1, . . . , n} exists such that if = j then 

1. either Sgn(ReXi) = Sgn(ReSj) 

or ReXi = ReSj — 0, but ImXi ^ and ImSi ^ 
or Xi = Sj = . 

2. the algebraic and geometric multiplicities of Xi are equal to those of 5j, 



for all i G {1, . . . , n}. 
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We will use the notation of section [Ol to denote an eigenvalue configuration. For example the eigenvalue 
configuration f}j^ denotes the equivalence class {(is\, —is\, —S2 + IS2, — S2 — IS2) | 8\, S2, S3 € R>o}- 

Our starting point here is the homogeneous reduced centralizer unfolding from section f4. 21 

Hr{v) = v x M x + V2M4 + v z hh + v 4 M$ + v b M b . 

We look for eigenvalue configurations: fjj, 0%02, 13(3, 1 . Also see tabled] However we must keep in 
mind that Ti. r is a reduced unfolding of L provided that v b 7^ 0. The following result characterizes the 
boundary of the stability domain and moreover shows that z/5 is unimportant when considering eigenvalue 
configurations only. 

Lemma 5.1. The eigenvalue configurations ofTL r {v) are of the types listed when the parameters v satisfy 
the conditions 



Configuration 


Conditions on v 


Pi 


Wi - vl)Wt + v't) + 


v x v i — 0? ^1 = ^4 = then v\ < | 


I3xl3 2 


vx = 0, V2VX = 




13(3, (3 2 


v x = 0, v± = 0, v 2 = 


±z/ 3 



For all cases we need that v b 7^ 0. 

Proof. In the proof of this lemma we will use the construction described in section f3. 31 Here the map tp 
parameterizes the hyper-surface of coefficients of polynomials having roots a±if3, ±ij. Such polynomials 
are of the form q(t) — ((t - a) 2 + l3 2 )(t 2 + j 2 ). Thus we have 

i/) : R 3 -+ R 5 : (a, f3, 7) ^ ((a 2 + /3 2 ) 7 2 , -2a 7 2 , a 2 + (3 2 + 7 2 , -2a, 1). 

Eliminating (a, (3, 7) from ij}(a,(3,j) — (00,01,02,03,04) we find /(a) = aoa§ + 040^ — 040203. Then for 
all t 6 R>o we have /(to) = t 3 f(a) and /(i 4 ao, i 3 ai, i 2 02, ta3, 04) = t 6 f(aQ, ax, 02, 03, 04). The map tj> 
maps the unfolding parameters v in Ti. r {v) to coefficients of the eigenvalue polynomial of 7i r (v), The 
explicit expression of (j>(v) is rather involved, therefore it will be omitted. Imposing the condition that 
H-riv) has eigenvalues of the form a ± i(3, ±i 7 we require that G{v) = where G = f o <fi. From this last 
equation we infer that G is homogeneous of degree 6. After some computations we find 

G(u) = f(4>(u)) = -64[^ 2 + vl] [{ul - v 2 2){v\ + v 2 ) + vtvt] . (8) 

The first factor in this expression is different from zero because ^5 is non-zero. This means that the 
boundary of the stability domain is in fact defined by F{v) = where F : M 4 -> R and 

F{v) = (y 2 - v 2 ) {v 2 + v 2 ) + v 2 v 2 . (9) 

Note that F is homogeneous of degree 4. A closer inspection reveals that we have to add the inequality 
v 2 I if v i — v i — 0- Moreover the set implicitly defined by F(v) — contains the boundary of the 
stability domain because we did not impose conditions on a in the eigenvalues a ± if}. H 

Remark 5.1. In [3] Bottema looks at singularities of the hyper surface aga 2 + a 2 — 040203 = in the 
space of coefficients of eigenvalue polynomials. This turned out to be sufficient to explain the so called 
destabilization paradox, see |29l 145] . ► 

Remark 5.2. Setting ao = 1 in aoo|+a 2 — 040203 = we get one of the Steiner surfaces: a§+af — 040203 
0. see [12 • This surface has two Whitney umbrella singularities. ► 

From now on we assume that v§ has a fixed non-zero value. With some abuse of notation we take i/£l 4 . 
The zero set of the function F in equation ^ defines a hyper surface. 

Definition 5.2. The set C = {y G R 4 | F(v) — 0} will be called the critical set. 

Since we only required that the eigenvalues of TL r (y) are a±i[3, ±17 without imposing the condition that 
a < 0, the critical set contains the boundary of the stability domain. In the next sections we proceed 
as follows. First we determine the local and global properties of C. Using the function F we find the 
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singularities of C. Because of the homogeneity of F they come in straight lines emanating from the origin 
(in K 4 ). This allows us to further reduce the problem to a 3-sphere transversally intersecting the critical 
set C in what will be called the critical surface S. The local properties of S can almost immediately 
be read off from those of C. For the global properties of S we use a Whitney stratification of S which 
we extend to a Whitney stratification of the 3-sphere. On the strata the eigenvalue configuration is 
constant. This allows us to identify the stability domain on the 3-sphere and describe the singularities 
on the boundary. More generally we thus obtain a description of a small neighborhood of L in gl(4,R). 



5.2 Local and global properties of the boundary of the stability domain 
Singularities of the critical set C 



We start with several properties of F in definition 15.21 Here we find the location of singular points of C. 
It is more convenient to characterize them as critical points of the critical surface to be defined shortly. 

Lemma 5.2. The polynomial F has the following properties 

a) F is homogeneous: F(tv) — t A F{v) for iel. 

b) For Vi e {-1, 1}, F{axv\, ■ ■ ■ ,04^4) = F{u x , . . .,i>i). 

c) Critical points of F in C come at least in straight lines because of the homogeneity. They are 
{(0,0, s,t) \s,te M} and {(0,s,t,0) \s,t€ K}. We furthermore find 



HessF(0, 0, s, t) 



(2{s 2 +t 2 ) 
-2t 2 

V 




0/ 



These points are transverse self-intersections except when t = 0. Later on we will relate t — to 
an intersection of self-intersections. For the other points we have 



/2(s z - - / 



HessF{0 7 s, t, 0) = 



2 _ +2\ 




V 





-2s 2 / 



Again we find transverse self-intersections when t 2 > s 2 . There are degeneracies at t — ±s that we 
will relate to Whitney umbrellas. 



Further reduction to critical surface 

To further simplify the analysis we use the fact that F is homogeneous. This strongly suggests to restrict 
to the intersection of C with the 3-sphere given by \\v\ \ = 1. Note that C is transverse to this 3-sphere 
so that we will not introduce intersection singularities. Thus we consider 

F{v) = {v 2 -vl){vi + vt)+v 2 vl=Q 

V 2 +tf +f| +z/f = 1, 

which is equivalent with 

\G{v)=vl-2v\vl-vlvl = K 

\ V 2 + v 2 + 4 + v 2 = l [ 1 

By eliminating ^3 from the equation F(v) = we obtain G(v) = v\ — 2v\v 2 — v\v\ = 0. We give 
the following geometric meaning to this manipulation. The 3-sphere ||^|| = 1 can be considered as two 
3-discs D ± = {(vi,V2,Vi) e M 3 \ v\ + v 2 + v\ < 1, 1^3 = ±\/l — (y\ + v\ + ^|)} glued smoothly along 
their common boundary = 0. 
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Definition 5.3. In the 3-sphere we call S = {v € R 4 | G(v) = 0, = 1} the critical surface. 

The next result immediately follows from the homogeneity of F . 
Corollary 5.3. The set C is a cone over S. 

Therefore we concentrate on the critical surface S. It consists of two surfaces S^, where ± refers to the 
sign of 1/3. Here are the 2-dimensional parts of {v S \ G(v) = 0}. Due to the symmetries of P, 
S~ is a copy of S + . The critical surface determines a decomposition of the 3-sphere in 3-dimensional 
open parts where we have different eigenvalue configurations. The latter is the topic of section [5] To 
obtain this decomposition we use the Whitney stratification of the critical surface and the corresponding 
incidence diagram. 

Singularities of the critical surface S 

The next proposition describes the local structure, in particular the singularities of S. The notation 
Pi, . . . , P4 and L\. .... L§ will become clear when we parameterize S + and S~ later on. 

Proposition 5.4. The critical surface has the following singularities 

a) Four lines of simple self-intersections: (yi, z/j, fi) — (0, 0, t) for t G [—1, 0) U (0, 1] and (yi, i>2, ^4) = 
(0,t,0) for t G [-■ l-s/2,0) U (0, i\/2], labelled L 1 ,L 2 ,L 5 ,L 6 for S + and L 3 ,L 4 ,L 5 ,L 6 for S~, 

b) Two points of self-tangencies at (^1,^2,^4) = (0,0,0) where the lines of self-intersection meet, 
labelled P5 and Pq for S + and S~ respectively, 

c) Four Whitney umbrella points at (fj., 1*2,1/4) = (0, ±\"/2, 0), labelled P\,Pi for S + and P 3 ,P4 for 

s-. 

Proof. The location of the singularities follows from lemma 15.21 Here we prove the nature of the 
singularities in parts [b] and First we consider G, see definition 15.31 locally at the point (0,0,0): 
G(x, y, z) = x 2 — 2x 2 y 2 — y 2 z 2 . The change of coordinates £ = x\J\ — 2y 2 , r] = y and ( = z yields that 
locally G(x, y,z) = is equivalent to £ 2 — rfC, 2 — (or even £(£ — nQ) = in yet another set of local 
coordinates). 

Let us next consider G locally at the point (0, |V2, 0) then a local change of coordinates yields the 
standard form of the Whitney umbrella. Indeed 

G(x, fV2 -y,z) = x 2 y(2V2 - 2y) - z 2 {\ - V2y + y 2 ) 

then by the change of coordinates 

i = 2x^l-\V2y, n=±V2y, C = z^ 2 -V2y + y 2 
equation G(x, \\[2 — y, z) — is equivalent to ^ 2 n — ( 2 — 0. ■ 

Global properties of the critical surface S 

We begin with the observation that the critical surface S as defined by equation (fT0| is a ruled surface. 
To see this, write the first part of equation (fTU]) as (1 — v\)v\ — v\v\ = and disregard the second part for 
the moment. Then S is formed by lines through the ^2-axis parallel to the v\, ^4-plane. In this respect S 
is a Catalan surface and in particular a right conoid. Again from the equation, but more easily from the 
parameterization below, we infer that S is equivalent to one of the Pliicker family of conoids, namely for 
n = l, where n is the index of the family, see [3l I18j. 

Imposing the additional condition that ||z/|| = 1, we find the following parameterization of S. 

Proposition 5.5. The critical surface S consists of two surfaces S being the images of [—1, 1] x [0, 2n] 
in D ± for the map (f> . This map is defined as 

± : [-1, 1] x [0,2tt] -> S ± : (s,t) i-> {\V2s cos t, \V2 cos t,s sin t). (11) 

It is one-to-one in most points, the exceptions are ^(s, = ± (— s, ^) for s 6 [— 1, 1] and (/) ± (0,<) = 
</> ± (0, 2n — t) for t G [0, 2tt], where ± is two-to-one. 
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Proof of proposition 15.51 It is easily checked that for all (s, t) e [—1,1] X [0, 2-k] we have G(</> ± (s, t)) = 0. 
Conversely, let v satisfy G(y) = 0. Fix V2 = c then G{v) = reduces to (1 — 2c 2 )z/ 2 — c 2 v\ = which 
clearly has pairs of straight lines as solutions as long as 1 — 2c 2 > and c > (a ruled surface). For 
fixed t, ± (s, t) maps onto these straight lines. If c = or 2c 2 = 1 the solutions are single straight lines 
which we find for t = \ir and s — respectively. Thus (fr^ maps onto the 2-dimensional part of S. The 
properties ± (s, f ) = ± (— s, ^L) and ^(0, t) = ^(O, 2tt — t) are again easily checked. ■ 

Note that the points where the map (jy^ fails to be one-to-one correspond to curves of self-intersection of 
S. See figure [3] for the domains of the maps ± . With help of the previous proposition we find the global 
structure of S. To this end we use a Whitney stratification of S. For sake of completeness we state a 
definition, see [2J. 

Definition 5.4. Let V as a smooth subset of a manifold M be a topological space. Furthermore let a 
collection of subspaces Ui C V, called strata, be given for i in an index set I . Then {Ui}i is called a 
Whitney stratification if 1) {Ui}i is a stratification and 2) for each pair X C Y and for all € Y with 
yk — > x € X , the tangent space T — lim^oo T yk Y satisfies T X X C T . 

From the differentiability properties of the maps <j> in equation (jlip we almost immediately get the 
following corollary. 

Corollary 5.6. The collection {Pi, ■ ■ ■ , Pq, L±, . . . , Lq, Si, ... , S^} forms a Whitney stratification of S . 
The organization of the stratification is shown in the incidence diagram in figure [^J 

There are several ways to describe the global structure of the critical set S. We start with the most 
'natural' one which leads us to a more useful description for identifying the stability domain on the 
3-sphere. 

Lemma 5.7. The open parts Si and S2 form a topological 2-sphere and similarly S3 and S4 form a 
topological 2-sphere. These two 2-spheres transversely intersect along the circle P5 U L5 U Pq U Lq except 
at the points P5 and P§. On each 2-sphere there are two (singular) crease intervals P\ U L\ U P5 U L2 U P2 
and P3 U L3 U Pq U L4 U P4. At these intervals the 2-spheres touch. The intersecting 2-spheres S\ U S% 
and S3 U S4 decompose the 3-sphere in the disjoint 3-dimensional open parts Vi, . . . , V4,. 

From the last property a more useful description in terms of Vi, . . . , V4 follows. The latter are domains 
in the 3-sphere with the same eigenvalue configuration. 

Lemma 5.8. The open parts S\ and Si form a topological 2-sphere S1US4 enclosing the open topological 
3-disc V\, similarly for the pairs (Si, 53), (S^S^) and (S^S^) enclosing V%, V3 and V4. 

Proof of corollary 15.61 and lemmas |5.7[ \ 5.8\ The incidence diagram can almost immediately be read off 
from the domains of </> ± in figure dJ For the sake of simplicity we use the same names for the strata in 
the domains as their images. Since is one-to-one in most points this should not lead to confusion. 
Points where ± is not one-to-one are indicated in figure [3] 

From the incidence diagram we infer that the points Pi, . . . , P4 and the lines L\, . . . , L4 are not essential 
for the global structure of S. Therefore the lines indicating their relation with other strata are dashed. 
Indeed, since Pi U L\ U P5 U L2 U P2 is a smooth interval in the domain of <fi + , so is its image. Thus 
<p + (Pi U L\ U P5 U L2 U P2) is a smooth curve with endpoints Pi and P2 (or to be formally correct </> + (Pi) 
and 4> + (P2)). This means we may shrink the lengths of L\ and L2 to zero without changing the global 
structure. Similarly for P 3 U L 3 U P 6 U L 4 U P 4 . 

Now we turn to Si U S^- From the domain of (f> + we see that S\ and S2 restricted to D + are smoothly 
attached along the lines L5 and Lq. Therefore they form a 2-disc in D + . The boundary of S\ and S% 
restricted to D + is a great circle, namely the image of the dashed line with four arrows in figure [31 This 
great circle is shown in figure [6] Similarly the restriction of Si and S% to D~ is also a 2-disc with the 
same boundary. These 2-discs are glued smoothly along their common boundary thus forming a 2-sphere. 
It is only a topological 2-sphere because there is a singular interval on each 2-disc. The same holds for 
S3 and S4. 

Looking again at the domains of (jy^ in figure [3] we infer that P5 U L5 U Pq U Lq is a smooth closed curve 
in S. In fact it is in S the boundary of each Si. Therefore Si U S4 for example is a topological 2-sphere. 
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The intersections of Si, S2, S3 and S4 with the boundary of D + , and of D~ , consist of two transversally 
intersecting great circles. Thus we see for example that S2 and S3 are both on the same side, let us agree 
that they are on the outside, of the topological 2-sphere Si U S4. Then the inside is an open 3-disc which 
we call V\. A similar construction holds for the pairs (Si, S3), (Ss,S2) and (S2,S4), yielding the open 
3-discs V2, V3 and V4. ■ 

6 Eigenvalue configurations on a neighborhood of L 

We are now in a position to describe the eigenvalue configurations on a small neighborhood of L in 
gl(4,R). In order to do so we apply a number of reductions. First we restrict to a neighborhood in the 
unfolding of L. By construction of the centralizer unfolding every matrix A near L is equivalent to some 
member of the centralizer unfolding of L, with p £ K 8 . Then, using the reduced unfolding, every 
element of C(p) is equivalent to some member of the reduced centralizer unfolding C r (v) of L, now with 
v £ M 5 . The third step consists of switching to the homogeneous reduced unfolding TL r (v). The results 
of the previous section show that for 7^ 0, a condition which is fulfilled in a small neighborhood of L, 
the eigenvalue configurations are constant along rays tv, with i/£l 4 and t £ (0, p]. Here p > depends 
on ^5. In view of remark [4.21 we may assume that ^5 is large enough so that we may take p = 1. Then 
we use homogeneity of eigenvalue configurations to restrict to the 3-sphere {v £ R 4 | \\v\\ = 1}. Thus we 
left with a three dimensional problem. 

For the latter we have a Whitney stratification {Pi, . . . , Pq, Li, . . . , L§, Si, ... , S4, Vj., . . . , V4} with in- 
creasing dimension, see figure [4] On each of these strata the eigenvalue configuration is constant. By 
choosing representative points and computing the eigenvalues we get table O From this table we infer 
that the stability domain of 7i r (y) is V3 . The boundary of the stability domain can be read off from the 
incidence diagram in figure [5J it consists of the components {S2, S3, Li, . . . , Lq, Pi, ... , Pq}. That is, the 
topological 2-discs S2 and S3 and moreover all singular lines and points, see figure El 

All that remains is to prove that table [5] is correct. 

Proof tabled is correct. The eigenvalue configurations follow from direct computation. The only thing 
left to prove is the difference between the nilpotent and the semi-simple case. The matrix L T (v) has a 
double pair of imaginary eigenvalues if v\ = v± = and ^3 = ±1^2- Then C r (0, V2, sv%, 0) — L + ^2(^4 + 
sMq) with s = ±1. By construction of the unfolding [L,M^ + sM^] — and using equation (O we see 
that (M 4 + sM 6 ) 2 = 0. Thus if S = L and N = z/ 2 (M 4 + sM 6 ) then S is semi-simple, is nilpotent, 
A 2 = and [S, A] = 0, so S + A is indeed the Jordan-Chevalley decomposition of £ r (0, v 2 , su^, 0). A 
double pair of imaginary eigenvalues only occurs at the points P\, . . . , P4 and C r is not semi-simple. I 
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